Scalelengths in Dark Matter Halos 
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ABSTRACT 

We investigate a hypothesis regarding the origin of the scalelengtli in halos formed in cosmo- 
logical N-body simulations. This hypothesis can be viewed as an extension of an earlier idea put 
forth by Merritt and Aguilar. Our findings suggest that a phenomenon related to the radial orbit 
instability is present in such halos and is responsible for density profile shapes. This instability 
sets a scalclcngth at which the velocity dispersion distribution changes rapidly from isotropic 
to radially anisotropic. This scalelength is reflected in the density distribution as the radius at 
which the density profile changes slope. We have tested the idea that radially dependent veloc- 
ity dispersion anisotropy leads to a break in density profile shape by manipulating the input of 
a semi-analytic model to imitate the velocity structure imposed by the radial orbit instability. 
Without such manipulation, halos formed are approximated by single power-law density profiles 
and isotropic velocity distributions. Halos formed with altered inputs display density distribu- 
tions featuring scalelengths and anisotropy profiles similar to those seen in cosmological N-body 
simulations. 

Subject headings: dark matter — galaxies:formation, evolution 



1. Introduction 

Some of the earliest studies to lay the foun- 
dation of galaxy formation have relied on ana- 
lytical methods to follow the coUisionless evolu- 
tion of dark matter {e.g., Gunn & Gott 1972; 
Gott 1975; Gunn 1977) and lead to pure power- 
law density distributions. Several other studies 
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have built on this early analytical framework. Fill- 
more & Goldreich (1983) and Bertschingcr (1984) 
have found similarity solutions for the spherical 
collapse of cosmological perturbations that pro- 
duce nearly pure power-law density profiles. The 
halos in these studies are populated solely by ra- 
dial orbits. Ryden & Gunn (1987) have intro- 
duced nonradial motion, but the resulting den- 
sity profiles arc still well-approximated by a single 
power-law over at least two orders of magnitude in 
radius, with noticeable deviations only on small 
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[log(r/rvir) < -2] and large [log(r/rvh.) > 0] 
scales (rvir is the virial radius). 

As computing power has increased, another ap- 
proach has been adopted to investigate cosmo- 
logical structure formation. N-body simulations. 
These studies do not suffer from the same limita- 
tions of adopted approximations that are inherent 
in the analytically-based work, however such simu- 
lations do rely on other approximations, like grav- 
itational softening. Given appropriate initial con- 
ditions, cosmological N-body simulations (CNS) 
can follow structure formation through non-linear 
development and investigate the impact of hierar- 
chical merging {e.g., Dubinski & Carlberg 1991; 
Navarro et al. 1996, 1997; Moore et al. 1998, 1999; 
Bullock et al. 2001; Power et al. 2003). These stud- 
ies agree that the equilibrium dark matter halos 
that form in CNS have nearly universal density 
profiles that are nothing like single power-laws. 
These halos generically have logarithmic slopes 
7 = —d\ogp/d\ogr that become larger with in- 
creasing radius; 7 « 1 at one- hundredth of a virial 
radius and 7^3 near the virial radius. There re- 
main, however, disagreements regarding the exact 
values and behavior of 7, especially at small radii. 

The fitting function that has been widely used 
to characterize the density profiles of CNS halos, 
the NFW profile (Navarro et al. 1996), has an ex- 
plicit scalelength that divides the inner and outer 
behaviors. The NFW profile 7 asymptotes to 1 as 
r — > and 3 as r 00. Recently, Navarro et al. 
(2004) have discussed another fitting function to 
describe the density profiles of higher resolution 
CNS halos. Unlike the NFW profile, 7 of the new 
function docs not approach asymptotic values at 
small and large radii and does not have the same 
well-defined scalelength as the NFW profile. In- 
stead, 7 of this new profile changes continuously 
with radius, implying that there are no regions 
in which the density profile behaves like a power- 
law. At the same time, another recent study (Die- 
mand et al. 2005) finds that the r ^ asymptotic 
power-law behavior persists in their high resolu- 
tion CNS. Given this disagreement, it is unclear 
whether or not CNS halos actually have preferred 
scalelengths. However, the NFW profile does pro- 
vide a reasonable description of density distribu- 
tions for galaxy and cluster scale CNS halos (Mer- 
ritt et al. 2005). For the sake of brevity, we will 
refer to the radius at which 7 = 2 as the scale- 



length. 

Despite these fairly minor disagreements re- 
garding the behavior of 7, CNS halos (as a class) 
appear qualitatively similar. This similarity does 
not extend to the analytically-based (AB) halos 
described above. Like previous studies that have 
sought to understand the differences between CNS 
and AB halos {e.g., Avila-Reesc, Firmani, & Her- 
nandez 1998; Lokas 2000; Nusser 2001), we are in- 
vestigating how the input physics differs between 
the methods and whether or not such physics can 
explain the differences. Wc find evidence that a 
phenomenon related to the radial orbit instabil- 
ity (ROI) can account for the differences as well 
as shed some light on the universality and scale- 
lengths of CNS halos. 

A brief introduction to the ROI is given in §2. 
We lay out the hypothesis, drawing links between 
the ROI, CNS halos, and AB halos, in §3. The 
demonstration of the hypothesis relies on evolving 
dark matter halos semi-analytically and is based 
on the method presented in Ryden & Gunn (1987). 
Our modifications to this method have been fully 
discussed in Williams et al. (2004). Briefiy, this 
method relies on spherical symmetry and ignores 
complications due to merging. This simplicity pro- 
vides us with a degree of control not present in 
CNS. At the same time, inclusion of secondary 
perturbations allows us to approximate the correct 
evolution of an isolated dark matter halo. Details 
of the method and results of testing the hypoth- 
esis are presented in §4. Wc discuss the link be- 
tween the ROI and the shapes of CNS halo density 
profiles in §5. The final section summarizes and 
presents the conclusions drawn from this work. 

2. Overview of the Radial Orbit Instabil- 

ity 

Early analytical and N-body studies have found 
that equilibrium spherical systems composed of 
purely radial orbits are unstable to forming bars 
or triaxial systems (Polyachenko 1981; Fridman & 
Polyachenko 1984; Binney & Tremaine 1987, §5.2). 
Although there is no general theory explaining the 
ROI, Palmer & Papaloizou (1987) argue that if a 
spherical system develops a small m = 2 distor- 
tion, precession of radial orbits can draw them to 
reinforce the distortion, thereby increasing the dis- 
tortion and leading to bar formation. Whatever 
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the exact mechanism is, it is clear that nonradial 
forces play a key role, as we will discuss in the next 
section. 

Mcrritt & Aguilar (1985) have further inves- 
tigated the ROI by creating equilibrium systems 
with initial anisotropy {(3=1 — cr^/cr) profiles 
described by, 

m = (1) 

where Va is referred to as the anisotropy radius 
(Merritt 1985). These profiles are completely 
isotropic (/3 = 0) for Tajr-Yji = cx), where 
is the initial half-mass radius. For smaller ra/ri/2 
values, the anisotropy increases from the center 
outwards. At the limit ra/ri/2 = 0, the en- 
tire model is composed of radial orbits. Mer- 
ritt & Aguilar (1985) list several fundamental re- 
sults of the ROI. First, an initially spherical sys- 
tem that undergoes the ROI transforms into a 
nearly prolate bar with long-to-short axis ratio 
sa 2 — 2.5. Second, for the anisotropy profile 
of Equation 1 there appears to be a fairly dis- 
tinct border between systems that form a sig- 
nificant bar (ra/ri/2 < 0.2) and those that re- 
main relatively spherical {ra/rii2 > 0.3); however, 
other anisotropy profiles produce different stabil- 
ity criteria. Third, the global anisotropy measure 
2Tr/Tt (where Tr and Tt are the kinetic energies 
associated with radial and tangential motions, re- 
spectively) is not a universal arbiter of whether or 
not a system is unstable; larger values of this ratio 
indicate instability, but there does not appear to 
be a definitive demarcation value (Polyachenko & 
Shukhman 1981; Barnes, Goodman, & Hut 1986). 
The Merritt & Aguilar (1985) work also highlights 
that the stability of any equilibrium system de- 
pends (at a minimum) on the anisotropy distribu- 
tion, with centrally isotropic systems being more 
stable than those with radially anisotropic cores. 
In brief, coUisionless systems will undergo the ROI 
if a sufficient fraction of orbits are predominantly 
radial. 

While we have so far discussed the ROI in equi- 
librium systems, coUisionless collapses also appear 
to relate to the ROI. The link with the ROI is that 
these cold collapses produce mostly radial motion, 
a condition favorable to the onset of the ROI. We 
compare the results of noncosmological and cos- 
mological collapses and discuss the signature of 



the ROI in the next section. 

3. Collapses &; the Radial Orbit Instability 

Studies of cold coUisionless collapses in non- 
cosmological settings (van Albada 1982; McGlynn 
1984; Aguilar & Merritt 1990; Trenti, Berlin, & 
van Albada 2005) report that the end results 
of collapses are good descriptions of observed 
surface brightness profiles of elliptical galaxies. 
For example, the collapses presented in van Al- 
bada (1982) result in mildly nonspherical systems 
with two-dimensional projected density profiles 
that are well-approximated by the Scrsic func- 
tion ln(S/I]e) = -bn[{R/ReY''" - 1] (Sersic 
1968), with n K, A like a de Vaucouleurs profile. 
These collapse products also display characteristic 
anisotropy profiles; the inner regions are isotropic 
and transitions to radial anisotropy occur at larger 
radii. Merritt & Aguilar (1985) also point out 
that such collapses are less centrally concentrated 
than purely radial collapses. We extend and fo- 
cus these thoughts by specifically investigating the 
link between density scalelengths and monotonic 
variations in velocity anisotropy. 

While some choices of initial conditions in the 
previously mentioned studies have an eye towards 
describing cosmological situations, their main 
goals relate to understanding elliptical galaxies. 
Two studies directly deal with the differences be- 
tween cosmological and noncosmological collapses. 
Katz (1991) argues that the ROI does not occur in 
cosmological settings because there is no correla- 
tion between final shape and initial 2T /W value. 
However, this study also reports that the final 
states of collapses have density profiles that are 
well-approximated by a Jaffe distribution (Jaffe 
1983), which is similar to a de Vaucouleurs (Sersic 
ri ~ 4) profile when projected to two dimen- 
sions. More importantly, it is stated that the 
equilibria are supported by anisotropic velocity 
dispersions, like the results of the noncosmologi- 
cal studies. Carpintero & Muzzio (1995) find that 
realistic cosmological collapses do not succumb to 
the bar instability, but, like Katz (1991), the fi- 
nal density profiles are approximated by a Jaffe 
distribution (again, similar to an n = 4 Sersic 
profile). Unfortunately, this study focuses solely 
on the final shapes of the collapses and does not 
present any kinematic information to allow a com- 
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parison of anisotropy profiles. Beyond these two 
studies, it is well-reported that halos formed in 
CNS tend to appear similar to the products of 
noncosmological collapse simulations; i.e., they 
appear prolate spheroidal/triaxial {e.g., Cole & 
Lacey 1996; Springcl, White, & Hernquist 2004) 
and have monotonically increasing anisotropy dis- 
tributions (Cole & Lacey 1996; Huss et al. 1999). 
All of these points lead us to conclude that cosmo- 
logical collapse simulations result in systems that 
have the same qualitative density and anisotropy 
profiles as those from noncosmological studies. 

The final states of these collapses are similar, 
but have the systems undergone the ROI? From 
the literature, it appears that most researchers 
consider bar formation to be the signpost of the 
instability. The Katz (1991) and Carpintcro & 
Muzzio (1995) studies suggest that the ROI is ab- 
sent in cosmological situations due to the lack of 
bars. Aguilar & Merritt (1990) refer to the merely 
slightly nonspherical results in van Albada (1982) 
as stable. 

We propose that bar formation may not be the 
only sign of the ROI, but possibly only the most 
flagrant. The end results of collapses like those in 
van Albada (1982), where no strong bar is formed, 
display a less drastic result of the ROI. The equi- 
librium states resulting from coUisionless collapses 
have similar anisotropy profiles; anisotropy in- 
creases (becomes more radial) with radius. This 
is in line with the discussion in Merritt & Aguilar 
(1985) which regards systems with isotropic cen- 
tral regions as more stable than those with radi- 
ally anisotropic central regions. From this point 
of view, the ROI changes the mostly radial initial 
velocity distribution into one that has an isotropic 
core surrounded by a radial "mantle" . As we con- 
tinue this discussion, we will use the term "mROI" 
(mild aspect of the ROI) to include the onset of 
the anisotropy profile and possible mild triaxial- 
ity and to distinguish it from the usual (and more 
extreme) bar-formation criterion. 

Figure 1 shows anisotropy profiles that can be 
constructed from information in several previous 
studies. Panel a) shows the results of van Albada 
(1982) (plus symbols) and the analytical expres- 
sion utihzed by Merritt & Aguilar (1985), Equa- 
tion 1. In the figure, we have set = Tg, where 
Te is the effective radius of the corresponding de 
Vaucoulcurs profile. Note that the analytical ex- 



pression is not a fit to the data points, it merely 
illustrates the initial anisotropy profile utilized in 
the ROI study of Merritt & Aguilar (1985). Panel 
b) shows the results from the "n = —2" (here, 
n refers to the power spectrum P{k) oc A:", not 
the Sersic profile) simulation of Cole & Lacey 
(1996) (asterisks) as well as the analytical expres- 
sion adopted by Carlberg et al. (1997), 



/3 = /3„ 



4r/ri78 



(2) 



where rns is the virial radius and we have cho- 
sen /?m = 0.3. This expression is intended to fit 
the data points from equilibrium systems. The 
last panel shows the results of the simulations of 
Huss et al. (1999). Admittedly, the amplitude of 
the anisotropy in the Cole & Lacey (1996) result 
(panel b) is much smaller than the other cases, 
but the monotonically increasing nature is still 
present. The mROI anisotropy profile is found 
generically in collapses. 

We find support for a link between the mROI 
and the presence of scalclcngths in CNS halos 
looking closely at the information from Huss et 
al. (1999). This work has investigated 5 models 
of isolated dark matter halo formation. Four of 
these have simple initial conditions that were cho- 
sen to highlight the impact of differing amounts 
of initial random motion. The fifth represents a 
standard NFW type halo. Their Model I is es- 
pecially interesting since it has no initial tangen- 
tial motions and tangential forces are not allowed, 
essentially forcing radial collapse. The resulting 
halo is described by a constant radial anisotropy 
{(3 = 1, diamonds in Figure Ic) and a single power- 
law density profile p (x r~^. Once this restriction 
on tangential motion is lifted (Models II-IV have 
increasing amounts of initial random motion and 
allow nonradial forces), the resulting anisotropy 
profiles show basically isotropic cores with (3 in- 
creasing with radius (Figure Ic), evidence that the 
mROI is present in these cases. The correspond- 
ing density profiles resemble NFW profiles, i.e., a 
scalelength is introduced. While generic equilib- 
ria do not have a unique relation between (3 and 
7, Hansen & Moore (2004) point out that a wide 
variety of CNS halos appear to have a specific link 
between the anisotropy distribution and 7. In §5, 
we further discuss this connection and the specific 
forms of the density and anisotropy profiles. 
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These previous studies provide evidence of 
two things. One, the mROL and its consequent 
anisotropy profile that has a distinct scalelength 
dividing an isotropic core from a radial mantle, is 
ubiquitous in collapses^. Two, this scalelength is 
linked to (if not the same as) the scalelength of 
the resulting density distributions. These points 
are the basis of the hypothesis that the mROI a) 
is the key physics missing from previous AB stud- 
ies of dark matter halo formation and b) leads to 
density profiles with scalelengths. We will now 
proceed to describe and discuss our testing of this 
hypothesis. 

4. Demonstrating the Link Between Anisotropy 
&: Density Profiles 

We are not interested in modeling the mROI. 
What we want to do is mimic the key result of 
the instability, the anisotropy profile, by suitable 
alterations of inputs. This approach will let us 
test the root of the hypothesis, the impact of the 
anisotropy distribution on the resulting density 
distribution. 

To create the halos we will investigate, we follow 
the procedure outlined in Ryden & Gunn (1987) 
and Williams et al. (2004). Initially thin spheri- 
cal shells expand with the Hubble flow until they 
reach their turn-around radius, at which time they 
collapse onto the mass interior. At the moment a 
shell begins to collapse, the integrated effects of 
secondary perturbations are introduced: a shell, 
which can be thought of as a mass of particles, 
is given a velocity dispersion. This can also be 
pictured as giving the shell a thickness, with apo- 
and pericentric distances. At the same time, the 
perturbation velocity is oriented randomly, intro- 
ducing angular momentum for that shell. The per- 
turbation velocities calculated per Ryden & Gunn 
(1987) are actually RMS values. Therefore, we 
follow Williams et al. (2004) and randomly choose 
the magnitude of the velocity from the Maxwell 
speed distribution (truncated at 4 standard devi- 
ations) with the RMS value from Ryden & Gunn 
(1987). Since our method relics on choosing ran- 

We note that this evidence cannot be considered conclusive 
as other studies have found that diflferont mechanisms, most 
notably density inhomogeneities, can also lead to similar 
outcomes (e.g., van Albada 1982; Londrillo, Messina, & 
StiaveUi 1991). 



dom magnitudes and directions, for the rest of this 
paper we will average over some number (usually 
20) of halos and refer to the average quantities as 
belonging to a single halo. 

The benefit of the WiUiams et al. (2004) for- 
malism is that it allows us to impose an arbitrary 
velocity dispersion profile and then study the halos 
that result from the subsequent collapse. However, 
before we explore variations in the anisotropy pro- 
file, we first establish a "standard" halo to use as 
a benchmark. By construction, it will be similar 
to the halos that have been created in Ryden & 
Gunn (1987) and WiUiams et al. (2004). The spa- 
tial density and phase-space density proxy p/a'^ 
profiles of the standard halo are shown in Figure 2 
(as usual, r2oo is the radius at which the halo den- 
sity is 200pcrit)- The top panel highlights the de- 
viations from p oc r~^. Over a substantial radial 
range (—2.5 < log(r/r2oo) ^ 0.0), the density is 
well described by a single power-law p cx r^^'®. 
This agrees well with the halo described in Ry- 
den & Gunn (1987, see their Figure 11). We show 
the p/ (7^ distribution for the standard halo in the 
bottom panel. This distribution is nicely approxi- 
mated by a power-law with exponent w — 2 over a 
wide radial range —2.5 < log(r/r2oo) < (Austin 
et al. 2005). We also present the velocity distribu- 
tions in Figure 3. The radial, one-dimensional tan- 
gential, and total velocity dispersion profiles along 
with the anisotropy distribution arc shown. Note 
that the velocities are basically isotropic over the 
entire range (the upturn outside r2oo is from in- 
falling material). Not surprisingly, given its scale- 
free power-law properties, we find that the stan- 
dard halo is a poor match to those from cosmologi- 
cal N-body simulations, represented in Figure 2 by 
a NFW profile of concentration 10 (dash-dotted 
line). 

We now turn our attention to altering the pa- 
rameters of the standard halo by varying the ve- 
locities caused by secondary perturbations. We 
introduce a factor v, which is a function of initial 
comoving radius x, to effect the desired changes. 
Since we are interested in mimicking the end re- 
sult of the mROI, wc introduce changes that pro- 
duce an anisotropy profile that is isotropic near 
the center and becomes radially anisotropic with 
increasing radius. 

The prescription for doing this in Williams et 
al. (2004) is a simple one; halve the perturba- 
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tion velocities in a specified comoving radial range 
(xio, Xhi) giving an inverted top-hat v distribution, 

{1 : X < xio 
0.5 : xio < X < Xhi (3) 
1 : X > Xhi. 

In this prescription, the perturbation velocities are 
simply multiplied by v. The impact of this pre- 
scription is to decrease the angular momentum 
imparted by perturbations in the specified radial 
range. This prescription also alters the total per- 
turbation velocity, thereby changing the energy 
imparted to the halo by the perturbations. The 
end result is a halo that loses the power-law prop- 
erties of the standard halo and instead resembles 
an NFW profile (see Figure 3 in Williams et al. 
2004). 

A representative v{x) distribution used in this 
study is shown in Figure 4. The distance measure 
of the i-axis is the initial comoving radius of a 
shell in Mpc. The algebraic expression is, 

v{x) = v*{x) - [v*{x) - d]e-(^-"^"'"/2'^', (4) 

where n and uj control the shape and width of the 
trough, d is the depth of the trough, Xh is the 
mid-point of the trough, 

v*{x) = G.5(i'o-t^oo)(l-tanh[2s(x-Xft)/(t^o-'^oo)])+i^( 

2s /{vq — Voo) controls the slope of v* at x/i, and 
Vq and Uao are the asymptotic values. These pa- 
rameters give us considerable flexibility to explore 
the impact of secondary perturbations on the final 
halos. After some experimentation, we have found 
that the interesting parameters to vary are vq, i^oo, 
and Xh- We fix the remaining parameters; n = 4, 
u> ~ 0.4, d = 0.5, s = 20. In general, i^o > 1 and 
Voo = 1- If we simply multiply the perturbation 
velocities by ly, then near the center {ly > 1) we in- 
troduce more angular momentum than in the stan- 
dard halo. Conversely, v < I saps angular momen- 
tum relative to the standard case. Through vari- 
ous trials, this change from ;/ > 1 to < 1 appears 
essential to reproducing the niROI anisotropy pro- 
file, and we point out that is simply a means of 
achieving this end. We have not determined the 
exact link between the behavior of ly, which is im- 
posed only at a shell's initial turn-around point, 
and the physical process of the mROI. It is possi- 
ble that the RMS perturbations calculated accord- 
ing to Ryden & Gunn (1987) are larger than those 



in CNS because of the one- dimensional nature of 
this method. If this is the case, then i/ < 1 is sim- 
ply reducing the effects of exaggerated secondary 
perturbations. While further work may elucidate 
the relationship between ly and the physics of the 
mROI, for the purposes of this study, it suffices 
that this V distribution gives rise to an niROI-type 
anisotropy profile. 

As it seems unlikely that the mROI changes the 
energies of particles, we turn to a new prescrip- 
tion that leaves shell energies unchanged and in- 
stead affects the orientations of the perturbation 
velocities. The prescription presented here (and 
shown in Figure 4) uses the following i> distribu- 
tion parameters; = 1.6, Voc = 1-0, x^ = 0.8 
Mpc. The magnitude of the perturbation veloc- 
ity remains unchanged from its value derived from 
the Maxwell speed distribution, but i/ > 1 in- 
creases the probability that the velocity is ori- 
ented tangentially and vice versa. We reiterate 
that this prescription is intended solely to produce 
an mROI-type anisotropy profile. There is no di- 
rect link between v and the physics of the mROI, 
and so the correspondence between the final pro- 
files is bound to be inexact. 

Figure 5 illustrates the resulting spatial den- 
sity and p/a^ profiles. We can see that adopting 
this new prescription has indeed provided a scale- 
length for the halo. While we have not set out 
to exactly recreate NFW halos, the similarity is 
clear. Some of the differences between the average 
and NFW profiles is due to the fact that many in- 
dividual halos are good approximations to NFW 
profiles with different concentrations (5 < c < 9). 
Averaging these together tends to smear the peaks 
into a broader profile. Looking at the distribution 
in panel (b), we see a profile that is a decent ap- 
proximation to a power-law over a couple orders 
of magnitude. It also appears to be a bit shallower 
(slope « —1.8) than the standard halo's p/c'^ pro- 
file and more in line with the results of Taylor 
& Navarro (2001) who find a slope of —1.875 for 
three CNS halos. Figure 6 contains the same ve- 
locity information as Figure 3. The break in the 
power-law that is present in the density has an 
analogous presence in the velocity distributions. 
Also, the anisotropy profile is familiar from the 
preceding discussion of the mROI; i.e., the veloci- 
ties are isotropic near the center with a fairly sharp 
transition to anisotropy. This result is in qualita- 
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tive agreement with the discussion of collapses in 
Merritt & Aguilar (1985); the ROI reduces the 
central density relative to purely radial collapses. 

To this point, the discussion and figures pre- 
sented in this section are all based on a single set 
of parameters that determine the v distribution. 
We have done some parameter space exploration 
and have found that the parameter can have 
an important impact on the final profiles. This is 
not very surprising, since it primarily determines 
the impact of the isotropic core. As vq is reduced 
from 1.6 to 1, the inner 7 of the density profile 
changes from « 1 to « 1.3. Also, while we are not 
trying to exactly reproduce NFW halos, we have 
noted that by eliminating the introduction of ran- 
dom perturbation velocity magnitudes (i.e., just 
using the RMS values) we can make the agree- 
ment between our halos (not shown) and NFW 
halos better. We speculate that this occurs be- 
cause it provides "colder" initial conditions (in the 
sense discussed regarding the ROI). We also find 
that the anisotropy profile is hardly affected if v^a 
is varied between 1.0 and d (the effect of changing 
Voo is to lower the v value from 1.0 to d = 0.5 over 
the range 1.5 < x < 2.5, see Figure 4). 

A more thorough investigation has been made 
regarding the importance of Xh, the mid-point of 
the I' trough. This parameter directly impacts 
the extent of the isotropic core present in the fi- 
nal anisotropy profile. In Figure 7, we show den- 
sity profiles (left column) and the corresponding 
anisotropy profiles for 7 values of Xh surrounding 
the canonical value; from top to bottom Xh—{0.6, 
0.7, 0.8, 0.9, 1.0, 1.1, 1.2}. We note that all 
the density profiles have scalclengths, but the 
agreement with the NFW form is within the er- 
ror bars only for Xh = {0.7,0.8,0.9}. More im- 
portantly, these models highlight the correspon- 
dence between the density scalelength and the 
anisotropy scalelength rj,, defined here to be the 
radius at which (3 = 0.5, as in Equation 1. Fig- 
ure 8 displays this correlation (solid line) and also 
shows how Ti, the largest radius at which (3 — 0, 
changes with (dash-dotted line). The density 
and anisotropy scalclengths are roughly the same 
(within « 20%) in these halos, further support- 
ing the hypothesis that mROI induced anisotropy 
profiles lead to scalelengths in density profiles. 



5. Relating the mROI and Shapes of N- 
body Halo Density Profiles 

In this section, we take up a discussion of the 
"true" shape of N-body halos. The preceding dis- 
cussion has purposely ignored the differences be- 
tween the "canonical" NFW profile and other CNS 
halo fitting functions {e.g., Navarro et al. 2004). 
One must realize that our lack of detailed un- 
derstanding of the physics of halo formation has 
forced us to approach halo structure from a very 
reactive posture; halos are formed and then func- 
tions arc chosen that resemble the resulting den- 
sity profiles. Unfortunately, these functions are 
somewhat arbitrary and do not, by themselves, 
lead to a deeper understanding of halo dynamics. 
The central idea of this study, that the mROI leads 
to an anisotropy distribution that is linked to the 
density distribution, can potentially provide a way 
around this arbitrariness. A full understanding of 
this physical mechanism can provide us with an 
expected density profile for collapse equilibria. In 
reality, the ROI is not fully understood and so no 
ab initio prediction can be made. At present, we 
are left to continue the passive role of finding the 
function(s) that best explain the simulation results 
and we now discuss the density profile specifically. 

The density profile that has long been found as- 
sociated with both noncosmological {e.g., van Al- 
bada 1982; Londrillo, Messina, & StiaveUi 1991) 
and cosmological {e.g., Katz 1991; Carpintero & 
Muzzio 1995) N-body simulations is the de Vau- 
couleurs (Sersic n ~ 4) profile. However, Sersic 
profiles are not very similar to the NFW profile, 
the dominant density fitting function over the past 
decade. While NFW profiles do fit CNS halo den- 
sity profiles well, recent work by Navarro et al. 
(2004) suggests that the asymptotic behavior of 
the NFW profile as r — > and r — » 00 does not 
fit density profiles as well as a function similar to 
the Sersic form (however, for a dissenting view see 
Diemand et al. 2005). Additionally, Dalcanton & 
Hogan (2001) and Merritt et al. (2005) have pre- 
sented evidence that the Sersic form is closer to 
the true shapes of halo density profiles. As the 
density profiles of the various ROI studies listed 
in this paper also find Sersic (specifically de Vau- 
couleurs) profiles, the N-body findings mesh nicely 
with the hypothesis that the ROI is a key part of 
forming collapse equilibria. 
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Taking the hnk between the niROI and cqui- 
hbrium density distribution as estabHshcd, we 
now speculate about the causal relationship be- 
tween them. Kazantzidis et al. (2004) point out 
that there are a wide variety of equilibria avail- 
able to halos with density profiles like those from 
CNS. This variety arises from different velocity 
distributions, but we have shown that CNS ha- 
los display mROI-type anisotropy profiles exclu- 
sively (see also Hansen & Moore 2004). Wc pro- 
pose that this choice between the equilibria is dy- 
namically set by the mROI, which brings about 
the anisotropy distribution that is associated with 
a Sersic density profile. As mentioned before, we 
cannot currently demonstrate this association ana- 
lytically, but Trenti, Bertin, & van Albada (2005), 
following the work of Stiavelli & Bertin (1987) and 
Merritt, Tremaine, & Johnstone (1989), discuss 
how the Sersic density profile goes hand-in-hand 
with the isotropic core and radially anisotropic 
mantle associated with the mROI. They present 
a distribution function meant to describe sys- 
tems that have experienced incomplete violent re- 
laxation. When this distribution function is in- 
tegrated over velocity space, the density profile 
is very similar to the Sersic form. Integrating 
this function over spatial coordinates, one derives 
an anisotropy profile that has the characteristic 
mROI shape. In the context of this idea, the link 
between anisotropy and density profiles is present 
regardless of initial conditions (i.e., cosmological 
or noncosmological), but we speculate that such 
different initial conditions most likely lead to dif- 
fering Sersic n values. Returning to the "true" 
shape of N-body halos, we conclude that Sersic, 
rather than NFW, profiles are better descriptions 
of the density profiles of CNS halos and that a 
fuller understanding of the ROI will provide in- 
sight to the actual density profiles. 

6. Summary &: Conclusions 

There is an important difference between mod- 
els of coUisionless dark matter halos created from 
analytically-based methods and those resulting 
from cosmological N-body simulations. Unless the 
input physics (Avila-Reese, Firmani, & Hernandez 
1998; Lokas 2000; Nusser 2001) or parameters are 
varied (as in §4), the halos that result from ana- 
lytical calculations tend to be best described by a 
single power-law density profile and approximately 



constant isotropic velocity dispersion anisotropy 
profiles. Cosmological N-body simulation halos 
tend to have density distributions that steepen 
and anisotropy profiles that become more radial 
with increasing radius. 

We have explored the hypothesis, akin to the 
discussion in Merritt & Aguilar (1985, §6), that 
the presence of scalclcngths in cosmological N- 
body simulations can be attributed to a mild as- 
pect of the radial orbit instability, which produces 
equilibrium halos with isotropic cores surrounded 
by regions of radial anisotropy. The hypothesis is 
that it is this fundamental change in the charac- 
ter of the orbits supporting the halo that leads to 
scalclcngths. Cosmological N-body simulation ha- 
los incorporate the physics of the mild aspect of 
the radial orbit instability and present such scale- 
lengths. Analytically-based halos do not have this 
instability "built-in" and therefore lack the conse- 
quent scalclcngths. 

Utilizing an extension of the analytically-based 
method of Ryden & Gunn (1987), wc have inves- 
tigated whether or not a preferred scalelength can 
be introduced into a halo that would naturally 
have a single power-law density distribution. To 
mimic the impact of the mild aspect of the radial 
orbit instability, we have artificially altered the ve- 
locities introduced by secondary perturbations to 
be more tangential near the center and more ra- 
dial in the outer parts of a halo. This is an ap- 
proximate way to include the effect of the instabil- 
ity in a technique that does not explicitly include 
such physics. In our prescription, the magnitudes 
of the perturbation velocities arc left unchanged, 
only the directions arc altered. This leaves the 
energy budget of the halo unchanged, which we 
argue is an acceptable approximation to the insta- 
bility. The important result is that the final ha- 
los have an anisotropy distribution reminiscent of 
those from cosmological N-body simulation halos 
and a density profile akin to an NFW profile. An 
exploration of the parameter space of our models 
shows that the presence of a scalelength is insensi- 
tive to the relative sizes of the isotropic and radi- 
ally anisotropic regions, and the anisotropy radii 
of halos are approximately equal to their density 
scalelengths. 

If cosmological N-body simulations are environ- 
ments in which the mild aspect of the radial orbit 
instability plays a significant role (and they ap- 
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pear to be so), two central features of such simula- 
tions have a ready explanation. One is the seeming 
universality of halo density profiles which appears 
in both noncosmological and cosmological N-body 
simulations as Sersic density profiles. A physical 
mechanism, like the mild aspect of the radial or- 
bit instability, that generically acts in collapses, 
even when mergers are absent, provides a simple 
generator for such universality. The outcome of 
the instability is a variable anisotropy distribu- 
tion with isotropic orbits in central regions and 
radially anisotropic orbits in outer regions. This 
anisotropy distribution gives rise to density pro- 
files with scalelengths, the second defining feature 
of halos formed in cosmological N-body simula- 
tions. 
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Fig. 1. — Plots showing the behavior of the 
anisotropy parameter (3 for various simulations, a) 
The plus symbols denote the results of van Albada 
(1982). The line is the initial anisotropy distri- 
bution utilized in Merritt & Aguilar (1985) (see 
Equation 1). b) The asterisks mark the values re- 
sulting in the n = — 2 simulation of Cole & Lacey 
(1996). This line is the distribution suggested in 
Carlberg et al. (1997). c) The following symbols 
represent the various models described in Huss et 
al. (1999); diamonds - Model I (purely radial), tri- 
angles - Model II, squares - Model III, crosses - 
Model IV, filled circles - Model V. 
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This 2-column preprint was prepared with the AAS lAT^X 
macros v5.2. 
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Fig. 2. — Standard halo spatial and phase-space 
density proxy distributions, a) Multiplying the 
spatial density by enhances the deviations from 
that power-law. The dash-dotted line is an NFW 
profile with a concentration of 10. b) The phase- 
space density proxy distribution. The \-a error 
bars are shown and the normalization is arbitrary 
in both panels. 
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Fig. 3. — Standard halo velocity distributions, a) 
The radial velocity dispersion profile, b) The one- 
dimensional tangential velocity dispersion profile, 
c) The total velocity dispersion profile, d) The 
anisotropy profile. The dash-dotted line shows the 
anisotropy of purely radial orbits. The solid line is 
for isotropic orbits. The l-cr error bars are shown 
in all panels. 



Fig. 5. — The spatial (a) and p/a^ (b) density 
distributions resulting from the prescription de- 
scribed in §4. The dash-dotted line in (a) is an 
NFW profile with a concentration of 8. The 1- 
(T error bars are shown and the normalization is 
arbitrary in both panels. 
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Fig. 4. — The distribution of v values versus initial 
comoving radius in Mpc. See §4 for the utility of 
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Fig. 6. — The velocity distributions resulting from 
the prescription described in §4. The panels are 
the same as in Figure 3. 
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Fig. 7. — Density (left column) and anisotropy 
(right column) profiles for different values of Xh- 
From top to bottom, x-/i={0.6, 0.7, 0.8, 0.9, 1.0, 
1.1, 1.2}. The dash-dotted line in the density pro- 
files is an NFW profile with a concentration of 
8. The ordinate values listed are the same for 
each panel in their respective columns. The solid 
and dash-dotted lines in the anisotropy profiles 
mark isotropy {(3 = 0) and total radial anisotropy 
(/3 = 1), respectively. Each model is the average 
of 10 individual halos. 




Fig. 8. — The correlations between anisotropy ra- 
dius rb (solid line) and maximum isotropic radius 
ri (dash-dotted line) and the density scalelength 
rd- The dashed line represents a one-to-one corre- 
spondence and the numbers are the corresponding 
Xh values. 
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